Magnetized Strange Quark Matter and Magnetized Strange Quark Stars 
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Strange quark matter could be found in the core of neutron stars or forming strange quark stars. 
As is well known, these astrophysical objects are endowed with strong magnetic fields which affect 
the microscopic properties of matter and modify the macroscopic properties of the system. In 
this paper we study the role of a strong magnetic field in the thermodynamical properties of a 
magnetized degenerate strange quark gas, taking into account /3-equilibrium and charge neutrality. 
Quarks and electrons interact with the magnetic field via their electric charges and anomalous 
magnetic moments. In contrast to the magnetic field value of 10 19 G, obtained when anomalous 
magnetic moments are not taken into account, we find the upper bound B < 8.6 x 10 17 G, for the 
stability of the system. A phase transition could be hidden for fields greater than this value. 



I. INTRODUCTION 

There are two possibilities for the occurrence of a phase 
transition between hadronic and strange quark matter 
(SQM) that are well known. The first could occur at 
very high temperatures and very low baryon density in 
the early Universe, and the second, as suggested by Bod- 
mer [l| , at densities of higher order than the nuclear den- 
sity no ~ 0.16 fm -3 . This phase transition would occur 
in the Universe, every time that a massive star explodes 
as a supernova, with its consequent remnant. If Bod- 
mer's conjecture is true, SQM could be succeeded in the 
inner core of neutron stars where strange quarks would be 
produced through the weak processes with a dynamical 
chemical equilibrium among the constituents. It is also 
possible that after a supernova explosion its core forms 
directly a strange quark star (SQS) 

The key property of SQM is that it has a binding en- 
ergy that could be lower than that of 56 Fe over a rather 
wide range of the QCD parameters [1, Thus, it is 
worthwhile to seek connections between SQM or SQS 
that could explain the observations of anomalous radia- 
tions from anomalous X-pulsars (AXPs) and low energy 
7-ray radiation from soft gamma ray repeaters (SGRs) 
p. In several SQM studies the essential conclusion is 
that a more compact matter would be the cause of these 
observations Q • On the other hand, astrophysical obser- 
vations point out that compact objects are endowed with 
strong magnetic fields which should play an important 
role in neutron stars or SQS. It is believed that magnetic 
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fields larger than 10 14 G are the central engine of their 
radiations. 

As is well known, a magnetic field modifies the mi- 
croscopic properties of matter with the corresponding 
macroscopic implication. However, the role of the mag- 
netic field in SQM has not been fully studied and un- 
derstood. Ref. 8] is a pioneer work in this field. The 
thermodynamical properties of quark matter in a strong 
magnetic field have also been studied in Refs. us- 
ing the MIT bag model modified by the inclusion of the 
magnetic field in the Lagrangian. In these works, it has 
been confirmed that there is an anisotropy of pressures 
due to the strong magnetic field [ill, an d that the 
MIT bag model can be used satisfactorily to study the 
magnetized quark matter. A first approach to consider 
the role of the anomalous magnetic moment (AMM) for 
quark matter has been followed in Ref. 

There are theoretical and experimental studies which 
indicate that quarks have an AMM [H, [3 [Tj| [H, G3 . 
A stringent bound on the quark AMM has been ob- 
tained [161 ] from high precision measurements at LEP, 
SLC and HERA. Thus, the contribution from the AMM 
of the quarks could be significant in SQM. On the other 
hand, for electrons, the effects of the AMM turn out to be 
small over the range of magnetic fields typically attained 
in neutron stars (10 15 — 10 19 G) [H[ and, therefore, can 
be safely neglected in our calculations. 

The scope of the present paper is to study the role 
of the AMM in the spectrum of particles and its rele- 
vance in the thermodynamic properties and stability of 
the SQM, considering /3-equilibrium and charge neutral- 
ity. We consider that the constituents of SQM interact 
with the magnetic field via their charges and their AMM. 
The effects of this interaction with the AMM could be 
seen in two forms. The first one is that quarks and elec- 
trons have the energy quantized in Landau levels due to 
their charge. The normal magnetic moment is included 
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in the spectrum through an integer moment. This quan- 
tization leads to a description of the polarization of the 
particles and to a softer equation of state (EOS). The 
second one is due to the inclusion of AMM terms. This 
causes a further splitting of the energy levels and the 
spectrum of the particles is no longer degenerate. But 
the most relevant implication of AMM in the spectrum 
of the particles comes from the ground state: it could be 
zero depending on the magnetic field strength. Thus, the 
AMM for individual particles establishes upper bounds 
on the strength of the magnetic field. We analyze the 
consequence of this fact in the macroscopic properties of 
the SQM in /3-equilibrium. 

For individual constituents, we find a saturation value 
for the field that align particles parallel or antiparallel to 
the magnetic field, depending on the AMM sign. This 
value of the magnetic field corresponds to a maximum of 
spin polarization and to the magnetization independent 
of B. This is the usual paramagnetic behavior. Beyond 
this saturation field, the spin polarization, magnetiza- 
tion as well as other thermodynamical quantities, become 
complex, which points towards a ferromagnetism phase 
transition. 

Since our treatment is based on non-interacting parti- 
cles, we cannot address the question of whether this lim- 
iting value corresponds or not to a phase transition. To 
clarify this issue a detailed study of the spin-spin coupling 
becomes necessary. Let us also remark that the SQM in 
/3-equilibrium adds new restrictions to the upper bound 
on the magnetic field. In this case, the polarization of 
SQM depends on the individual polarization of the con- 
stituents and its orientation is related to the AMM sign. 
The system is so complex that the total polarization is 
not reached in any direction. Nevertheless, it is possible 
to find a critical value of the magnetic field beyond which 
the polarization becomes complex, lacking physical sense. 

Our results improve earlier works in three aspects. 
Firstly, we take into account Pauli paramagnetism in its 
relativistic version, because the one-particle energy is the 
solution of the Dirac equation including the AMM. This 
gives an important contribution to the physics of the sys- 
tem. The upper bound on the magnetic field for each par- 
ticle is lower than the one obtained classically. Secondly, 
we consider the anisotropics of the pressures within the 
MIT bag model. This leads to changes in the behavior of 
the total energy of the system for strong magnetic fields. 
Finally, all quarks are assumed to interact with the mag- 
netic field. The most important astrophysical implication 
of our study is the existence of a limiting value for the 
magnetic field. For SQM with electrons we find an upper 
bound on the magnetic field around 8.6 x 10 17 G. This 
allows us to conclude that there would not be quark stars 
with magnetic fields greater than this value. 

Another effect associated with the inclusion of the 
AMM in the particle spectrum is that it can avoid the di- 
vergence arising through the lowest Landau level (n = 0) 
in the calculation of the surface tension and curvature for 
dense matter, and which play a crucial role in the quark 



droplet nucleation process [8[. This means that it is still 
possible that a first-order phase transition to quark mat- 
ter can occur [l9| . A detailed study of this effect and 
its astrophysical consequences is in progress. The AMM 
might play a role if Bose-Einstein condensation due to 
the bosonization of fermions takes place. We expect in 
that case a ferromagnetic behavior, able to maintain the 
applied magnetic field self-consistently [2fj| . 

The paper is organized as follows. In Sec. [IT] we study 
the spectrum of constituents of the magnetized SQM: 
electrons and quarks. Sec. IIIII is devoted to the analy- 
sis of the thermodynamic properties of the magnetized 
SQM. In Sec. IIVI we establish the requirement for the 
stability of SQM in /3-equilibrium, and study the spin 
polarization and its implication for the thermodynami- 
cal properties of the system. In Sec. |V] we present our 
numerical results, including a comparison of the behavior 
of SQM when AMM are taken into account with the case 
when the latter are neglected. Finally, our conclusions 
are presented in Sec. EB 



II. 



SPECTRUM FOR THE CONSTITUENTS OF 
MAGNETIZED SQM 



The relativistic spectrum of electrons in the presence 
of a magnetic field with the inclusion of AMM is obtained 
from the Dirac Pauli equation 
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e is the electron charge, m e is the electron mass, F M „ 
is the electromagnetic tensor and a is the fine structure 
constant. Considering a constant uniform magnetic field 
B in the cc3-direction, the energy spectrum is given by 1 
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?/ = ±1 are the eigenvalues corresponding to the ori- 
entations of the particle magnetic moment, parallel or 
antiparallel to the magnetic field. 



1 Hereafter we use the units h - 
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An analogous equation can be written for quarks, so 
that the spectra for all the constituents of SQM have the 
form 
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i = (e,u,d, s), &i and rrii denote the charges and the 
masses of the particles, respectively. The quantities Qi 
are the corresponding AMM of the particles, 



Q e = 0.00116/x S: 
Q d = -0.97^ N , 



Q u = 1.85fi N . 

Q s = -0.58/xat . 
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where 



Mb = 
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2m e 
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2m T : 



5.79 x 10~ 15 MeV/G , 
: 3.15 x l(T 18 MeV/G. 



(8) 



In our calculations we shall take m u — rrid = 5 MeV 
and m s = 150 MeV for the light quark masses. The 
magnitudes of the so-called critical fields Bf (when par- 
ticle's cyclotron energy is comparable to its rest mass) are 
B% = 4.4 x 10 13 G, B c u = 6.3 x 10 16 G, B c d = 1.3 x 10 16 G 
and B c s = 1.1 x 10 19 G. 

It can be seen from the spectra (0 that, besides of the 
quantization of their orbits in the plane perpendicular to 
the magnetic field, charged particles with AMM undergo 
the splitting of the energy levels with the corresponding 
disappearance of the spectrum degeneracy. For the non- 
anomalous case, Qi = 0, the minimum energy is inde- 
pendent of the magnetic field strength and the magnetic 
field only quantizes the kinetic energy perpendicular to 
the field. In this situation, the energy is degenerate for 
Landau levels higher than zero. States with spin parallel 
or antiparallcl to the magnetic field orientation (77 = ±1) 
have the same energy. However, the anomalous case, 
Qi ^ 0, removes this degeneracy. In the latter case, the 
rest energy of the particles depends on the magnetic field 
strength. The ground state energy is 



Ei,o = mi (1 - yiB) . 



(9) 



The above equation leads to the appearance of a 
threshold value for the magnetic field at which the ef- 
fective mass vanishes, mi ~ |Qi|-B. The thresholds of the 
field, Bf = 1/yi, for all the constituents of the SQM are 
given by 

B s e = 7.6 x 10 16 G, B s u = 8.6 x 10 17 G, 

B s d = 1.6 X 10 18 G, B^ = 8.2xl0 19 G, (10) 



that are smaller than the ones obtained when the classical 
AMM contribution is considered [|[. 

The meaning of the ground state energy value for QED 
was discussed long time ago [2l|. Nevertheless, it was not 
emphasized enough the fact that, due to the degeneracy 
of the orbit center, such ground state level may be popu- 
lated by a larger number of particles. The expression © 
suggests that the energy of the particles becomes smaller 
than the corresponding one for the antiparticles, with 
the consequent creation of pairs. The sign of the energy 
state is an invariant property for particles and antiparti- 
cles. This also means that positive and negative energy 
levels of electrons will never cross each other, i.e. it exists 
a non-crossing property. The spontaneous pair creation 
in a magnetic field is forbidden. Thus, for individual par- 
ticles, the correct meaning of this "critical" field is that 
it corresponds to an upper bound. 

Let us remark that in the SQM scenario all the con- 
stituents interact with the magnetic field and are obliged 
to satisfy the equilibrium conditions. Under such con- 
straints, it turns out that the dominant threshold field 
comes from u quarks, thus leading to the upper bound 
B < 8.6 x 10 17 G (see section IVl below) . This result has 
an important astrophysical consequence, since the bound 
for SQM can be also extrapolated to the SQS scenario. 
If SQS exist, the maximum magnetic field strength that 
they could support would be around the above bound, 
i.e. 10 18 G. 



III. THERMODYNAMIC PROPERTIES OF 
MAGNETIZED SQM WITH AMM 

The MIT bag model is appropriate for the study of 
magnetized quark matter |9j. In that model, confine- 
ment is guaranteed by the bag and quarks are considered 
as a Fermi gas of noninteracting particles. Under these 
assumptions, it is possible to study the thermodynamical 
properties of a quark gas in a strong magnetic field. In 
this section we investigate the thermodynamical proper- 
ties of the SQM when the AMM is included. 

The inclusion of AMM and the consequent loss of de- 
generacy implies that the sum over Landau levels is re- 
placed by two sums 



(ii) 



71=0 



n—0 rj 



For each thermodynamical quantity, the summation over 
the spin orientation leads to two contributions, corre- 
sponding to particles with the spin aligned parallel or 
antiparallel to the magnetic field. Moreover, since parti- 
cles have positive or negative AMM, they have different 
preferences in the spin orientation with respect to the 
magnetic field. As we shall see below, this has important 
consequences in the EOS of the system. The most rele- 
vant comes from the lowest energy ground state, which 
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depends on the strength of the magnetic field and it could 
be zero (cf. Eq. ©). 

For a degenerate magnetized SQM, where only par- 
ticles contribute to the thermodynamical potential and 
temperature can be formally taken as zero, the expres- 
sion for the thermodynamical potential can be written in 
the form [j| 



fif = - Xi gf + hf 2 



diCiTTl^ 



4ir 2 



(12) 



where i = e,u,d,s and di is a degeneration parameter 
(d e = 2,d u ,d,s — 6). We have defined Xi = fj,i/m,i, where 
fii is the chemical potential; hi and gi are dimensionless 
functions given by 2 
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The sum over the Landau levels n is up to n % max given by 
the expression 



n % = T 

max 



(Xi + rjy.B) 2 - 1 



2B/B? 



(15) 



where I[z] denotes the integer part of z. 

The density of particles, defined as N = J2i^i with 



Ni = ff gives 
B 



2tt 2 
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phase transition takes place or not. The resulting ex- 
pression is 

M? 



B 



- +X t ), 
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where 
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(20) 



We notice that the magnetization and magnetic sus- 
ceptibility are not linear functions of B. Moreover, the 
requirement for the ground state energy Ei£ > is equiv- 
alent to the condition h^(n = 0) > (see Eq. (fT4"|)). Since 
all the thermodynamical quantities depend on h 1 ?, they 
loose their physical meaning when < 0. Therefore, 
the condition > means that the system cannot ad- 
mit a value of the magnetic field greater than B? , which 
indicates that some phase transition occurs. The mag- 
netization reaches a value independent of the magnetic 
field, as it occurs for a paramagnetic system, but beyond 
this limit it becomes complex and could be associated to 
a ferromagnetic transition. 



A. Pressure and energy density 

Let us now write down the expression for the 
anisotropy of pressures and for the energy density of SQM 
when the AMM is included. The energy density, U, for 
the gas of i-particles can be obtained from the energy- 
momentum tensor [9[ , leading to 



iNi 



(21) 



The magnetization of the SQM is M — J2iMi, with 
Mi = -f|. We find 



Mi = M° i Y l (Mt+Mi), 



Mf = 9 f Xi - (hf 2 + 2A±Tf ) ln 



with 



(2n + 1 - rf)B 



2B?^B/B?(2n+l-r]) + l 



- VViB ■ 



(17) 
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We may also calculate the magnetic susceptibility x — 
J2iXii which is defined as and can measure if a 



2 To simplify the notation, from now on we omit the Landau level 
subscript n in all quantities. 



Evaluating this expression we find 

u i = M i B y £,(ut + v; 



U±= Xq gt + ht 2 \n X -^ 



The pressures are obtained from the expressions 
P±=J2 P ^> Pi± = -CU-MiB, 

i 
i 

Using Eqs. (JT2J) and ^TTJ) we find 

P<\\=M i B y £(p+ l +^ , 

n 

P ± - x a ± h ± 2 ln Xl+ ^ 



(22) 



(23) 



(24) 



5 



and 



Pix=M^Bj2(Ptx + Pr±) > 
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IV. THE STABILITY CONDITION FOR 
MAGNETIZED SQM 

In this section we study the stability condition of the 
SQM in a strong magnetic field. In the context of the 
MIT bag model and in the absence of a magnetic field, the 
stability condition for SQM means to study the equation 



Pt + B, 



bag 



together with the total energy 



(26) 



(27) 



under the condition Pt = 0. 

In the presence of a strong magnetic field, the bag pres- 
sure Bb a g has an anisotropic form that depends on the 
B-direction in space, 



B 



bag 



B bag = y^.Pi -L 



(28) 



Since the magnetization is always a positive function, this 
anisotropy in the pressures implies Pj_ < Pu . Thus, the 
stability condition for strong fields changes from Pt = 
to P^r = or, equivalently, 



J2^ = -J2 M ^ B - 

i i 

The total energy becomes 

U T = ^(-MiB + x t Ni) . 



(29) 



(30) 



For weak fields P± = Pu and we get the expression 

U T = x i N i > ( 31 ) 

i 

which is in agreement with [22| where anisotropics due 
to strong magnetic fields have not been considered. 

A. Paramagnetism response: spin polarization 

In Ref. [8j the Landau diamagnetism related to charged 
particles in a magnetic field was studied, treating classi- 
cally the relativistic behavior and Pauli paramagnetism 



associated with the inclusion of AMM. Our starting view- 
point is different. We consider a relativistic equation of 
motion, taking into account both contributions from a 
relativistic point of view. In this sense, our treatment is 
more robust since the AMM is included in the spectrum 
of particles as a relativistic effect. 

As we have already shown, the thermodynamical quan- 
tities for each constituent of SQM have two terms, related 
to particles with spin up and down orientations. The 
number density is not excluded from this fact, so it is 
important to study its behavior because it gives us infor- 
mation about the spin polarization of the system. Let us 
rewrite Eq. (|16[) in the form 



Nt = Nj + N} . 



(32) 



In the absence of a magnetic field or when B — > we can 
see that Nj = N^. On the other hand, in the presence of 
a magnetic field, the relation for the number density im- 
plies the existence of a magnetic field strength threshold 
for which complete saturation of each constituent of the 
SQM occurs. Whether or not a complete saturation for 
all the particles involved in the system is attained will 
depend on the SQM equilibrium conditions. 
We can define the spin polarization rate as 
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(33) 



When Bf is reached, we have the following condition for 
each constituent particle, 



S 



l . 



s 



0. 



(34) 



which means that is maximum and = 0. The 
threshold field value Bf saturates the system and aligns 
all particles parallel or antiparallel to the magnetic field. 
This alignment depends on the sign of the AMM of each 
particle. From the thermodynamical point of view, this 
behavior could be understood as a paramagnetic response 
of the system and it makes an important difference. Be- 
yond this magnetic field value all the thermodynamical 
quantities become complex for a pure gas of particles. 

Let us recall that the conditions of /3-equilibrium and 
charge neutrality add new restrictions to the thresh- 
old values Bf. These values were computed numeri- 
cally for all the chemical potentials of the SQM con- 
stituents to study, in the next section, the spin polariza- 
tion in the regime of a strong magnetic field. As it turns 
out, the /^-equilibrium and charge neutrality are satis- 
fied only for magnetic field values below certain thresh- 
old. Above this, the number density of electrons and 
u quarks are fixed, independently of the magnetic field, 
and /3-equilibrium will require /i e < 0. Our numerical re- 
sults confirm that the main contribution to the threshold 
field comes from u quarks, which imply the upper bound 
B < Bf, = 8.6 x 10 17 G. 
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FIG. 1: (Color online) Chemical potentials fa for SQM as 
functions of the magnetic field strength B with (solid lines) 
and without (dashed lines) AMM. The vertical dot-dashed 
line corresponds to the threshold value ~ 8.6 x 10 17 G. 



FIG. 2: (Color online) The variation of the number densities 
with the magnetic field. The solid (dashed) lines correspond 
to the case with (without) AMM. The baryon density is fixed 
at the value ub = 0.4 fm -3 . 



V. SQM IN /3-EQUILIBRIUM AND CHARGE 
NEUTRALITY: NUMERICAL RESULTS 

In this section we perform a complete numerical study 
with the aim to determine all the relevant thermodynam- 
ical quantities for SQM and discuss its stability, taking 
into account /3-equilibrium and charge neutrality. This 
requires the solution of a system of equations to obtain 
the chemical potentials of all the species involved in the 
system. If SQM exists in the core of neutron stars or 
forms itself a SQS, weak processes will be responsible for 
the appearance of the s quarks. Once this occurs, the 
equilibrium among the constituents will be dynamically 
established. 

The three ingredients for the SQM in equilibrium are 
/3-equilibrium, charge neutrality and the conservation of 
the baryonic density ns- 

fi u + fl e = fid, M<i = Ms > 
2N U -N d -N s -3N e = 0, 

^(N u +N d + N s )=n B . (35) 

Here we assume that there is no neutrino trapping in 
the system, so that they do not play any role on the 
/3-equilibrium conditions. For a given baryon density 
(we take ns — 2.5 no = 0.4 fm~ 3 ) and magnetic field 
strength, these equations together with Eqs. (fTBj) allow 
us to determine the chemical potentials and evaluate all 
the thermodynamic properties of the system. At the end 
of this section we shall comment on the dependence of 
our results with the variation of ub ■ 

In Fig. [1] we show the chemical potentials Uj, i.e. the 
solution of Eqs. (|3"S")) , as functions of the magnetic field B. 
Note that the chemical potentials remain practically con- 
stant up to the threshold B s u around 8.6 x 10 17 G, which 
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FIG. 3: (Color online) Spin polarization for SQM as a function 
of the magnetic field strength B with (solid lines) and without 
(dashed lines) AMM. 



corresponds to the upper bound on the magnetic field 
determined by the /3-equilibrium condition. The varia- 
tion of the number densities for all the SQM constituents 
with the magnetic field B is shown in Fig. [2] Compar- 
ing this variation with and without the AMM inclusion, 
we can see that for relatively small values of the mag- 
netic field, B < 10 16 G, all the number densities remain 
practically constant. At around the magnetic field value 
of 4 x 10 16 G, the electron density A^ e start to increase 
with B, whereas the quark densities remain almost con- 
stant up to field strengths of 10 19 G. Above this value 
the s-quark density N s decreases with B and becomes 
negligibly small. The oscillations due to the presence of 
Landau levels can be seen for the case when the AMM 
is not considered. Clearly, the AMM curves are bounded 
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FIG. 4: (Color online) Magnetization as a function of the 
magnetic field with (solid line) and without (dashed line) 
AMM. 



FIG. 6: (Color online) Anisotropy of the SQM pressures as 
functions of the magnetic field strength. The cases with (solid 
lines) and without (dashed lines) AMM have been considered. 
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FIG. 5: (Color online) Susceptibility as a function of magnetic 
field for the two cases considered: without the inclusion of 
AMM (dashed line) and with their inclusion (solid line). 



FIG. 7: (Color online) Energy per baryon versus B. As in 
previous figures, the dashed line correspond to the case when 
the AMM is not included, whereas the solid line takes into 
account the AMM of all the SQM constituents. 



by the upper bound B^; for fields greater than this value 
all the number densities become complex. 

The total spin polarization of SQM, = J2t S p l > 
is plotted in Fig. [3] as a function of the magnetic field. 
We see that the spin polarization of the system increases 
with the increasing of the magnetic field. When the 
AMM are not included, a total polarization is achieved 
for B ~ 2 x 10 19 G. On the other hand, when the AMM 
are taken into account the system cannot reach a total 
spin polarization since for values greater than B^ the 
density number becomes complex. 

The behavior of the magnetization A4 is depicted in 
Fig. |4j It is always a positive quantity for fields greater 
than 10 16 G. It also exhibits the so-called de Haas-van 
Alphen oscillations, with increasing amplitude as B in- 
creases. This is even more noticeable when the AMM is 



not considered as higher values of B are allowed. The 
corresponding magnetic susceptibility \ is presented in 
Fig. [5] It shows the paramagnetism behavior of SQM 
for fields larger than 10 16 G. Below this magnetic field 
strength, \ has an oscillating behavior. The upper bound 
on the magnetic field around 8.6 x 10 17 G encloses a phase 
transition of second type, because at this value all par- 
ticles align parallel or antiparallel to the magnetic field 
with a positive value of the magnetization. 

Let us now consider the pressure and total energy of 
the system. We plot in Fig.[S]the behavior of the pressure 
with the variation of the magnetic field. The system 
persists being anisotropic when the AMM are considered. 
For vanishing AMM, the perpendicular component of the 
pressure P± goes to zero at about 2 x 10 19 G, when the 
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FIG. 8: (Color online) Energy per baryon versus ns with 
AMM included for B = 10 13 , 10 15 , 10 17 G. 
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FIG. 9: (Color online) Total spin polarization versus ns with 
AMM included for B = 10 13 , 10 15 , 10 17 G. 



total spin polarization is reached. On the other hand, the 
inclusion of the AMM forbids fields above the threshold 
Thus, for SQM the anisotropy in the pressures is 
relatively small, i.e, P± ~ Pn . 

Fig. [7] shows the behavior of the total energy per 
baryon with the magnetic field. We have plotted two 
curves: particles with AMM and without it. The figure 
confirms that SQM is stable up to the corresponding field 
threshold. For B < B^ the energy per baryon remains 
almost constant, and decreases for higher values of the 
magnetic field. 

To conclude this section, let us comment on the vari- 
ation of the thermodynamic properties with the baryon 
density Ub- In Fig. [5] we show the dependence of the total 
energy per baryon with the baryon density for three dif- 
ferent values of the magnetic field, 10 13 , 10 15 and 10 17 G, 
taking into account the AMM. As the baryon density 
increases, the total energy of the system increases too. 



It also turns out that SQM is energetically more sta- 
ble in the presence of a strong magnetic field 3 . For a 
fixed value of ns , we remark that there are no significant 
changes as the magnetic field varies (the three curves are 
almost indistinguishable, as can be seen from the figure). 
The same is observed for other thermodynamical quan- 
tities. Nevertheless, the total spin polarization increases 
with the increasing of the magnetic field strength and the 
decreasing of the baryon density ns, as can be seen in 

Fig.m 



VI. CONCLUSIONS 

We have studied the magnetized SQM in /3-equilibrium 
in the presence of a strong magnetic field. We have taken 
into account the Landau diamagnetism related to the 
quantization of the Landau levels as well as the Pauli 
paramagnetism, due to the presence of AMM for all the 
constituents of SQM. The influence of the paramagnetism 
in the system is more relevant than the diamagnetism be- 
cause it is responsible for the upper bound on the mag- 
netic field found for the SQM system. This bound is 
lower than the one obtained classically [8] . Furthermore, 
it implies that a phase transition should occur at this 
value, because all individual particles are aligned paral- 
lel or antiparallel to the magnetic field (depending on the 
AMM sign) in the ground state of the energy. 

For SQM in /^-equilibrium and with neutral charge the 
situation is mathematically complex. The condition of 
/3-equilibrium implies an upper bound on the magnetic 
field, B < B s u — 8.6 x 10 17 G. Above this value, the chem- 
ical potential of electrons becomes negative and all the 
thermodynamical quantities loose their physical mean- 
ing. As a consequence, a total spin polarization is not 
achieved, in contrast with the case without AMM, where 
such a polarization is reached for fields ~ 10 19 G. 

From the quantum statistical point of view the low- 
est energy states with AMM contain important physical 
consequences: for particles with mass rrii and anoma- 
lous magnetic moment Qi, the magnetic field has a crit- 
ical value given by the expression B? ~ mi/\Qi\. It re- 
mains to clarify if, for a magnetic field strength of this 
order, quark matter undergoes a phase transition. This 
question deserves particular study. On the other hand, 
it becomes clear that the stability condition of SQM is 
modified in the presence of a strong magnetic field. 

In this work we have shown the differences that the 
AMM introduces in all thermodynamic properties. We 
have concluded that magnetized SQM with AMM is sta- 
ble as it is in the case when no AMM is introduced. In 
both cases, magnetized SQM is more stable than SQM 



3 In the absence of a magnetic field, and for a given value of tib , the 
SQM energy per baryon is always higher than the value obtained 
when B ^ ||. 
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without a magnetic field. The pressures preserve the 
anisotropies found in [l2j for pure neutron and electron 
gases in strong magnetic fields. Nevertheless, a signifi- 
cant anisotropy cannot be reached due to the presence of 
the AMM. 

In summary, we have obtained a threshold value for 
the magnetic field which is equal to the saturation field 
for the u quarks. This value of the magnetic field is 
due to the restrictions of /3-equilibrium and it limits all 
the thermodynamical quantities. From the astrophysical 
point of view, our conclusions imply that, if SQS exist, 
they cannot support magnetic fields greater than 10 18 G. 

As mentioned before, if there is bosonization, which is 
otherwise expected, for instance, in the form of di- quarks, 
the model of Bose condensation developed in [2(1 [HI 
could be applied. For SQM, in that case, the ferromag- 
netic phase transition due to AMM would be guaranteed 
for fields of the order ~ 10 18 G. This could indicate an- 
other type of phase transition for SQM, in addition to 



the CFL [3 (or mCFL [H[) phases. 
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